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Abstract. We present medium resolution (R~1500) ISO-SWS 2.4-45/^m spectra of a sample of 29 galaxies 
with active nuclei. This data set is rich in fine structure emission lines tracing the narrow line regions and 
(circum-)nuclear star formation regions, and it provides a coherent spectroscopic reference for future extragalactic 
studies in the mid-infrared. We use the data set to briefly discuss the physical conditions in the narrow line regions 
(density, temperature, excitation, line profiles) and to test for possible differences between AGN sub-types. Our 
main focus is on new tools for determining the properties of dusty galaxies and on the AGN-starburst connection. 
We present mid-IR line ratio diagrams which can be used to identify composite (starburst + AGN) sources and 
to distinguish between emission excited by active nuclei and emission from (circum-nuclear) star forming regions. 
For instance, line ratios of high to low excitation lines like [O IV]25.9^tm/[Ne II]12.8/im, that have been used to 
probe for AGNs in dusty objects, can be examined in more detail and with better statistics now. In addition, we 
present two-dimensional diagnostic diagrams that are fully analogous to classical optical diagnostic diagrams, but 
better suited for objects with high extinction. Finally, we discuss correlations of mid-infrared line fluxes to the 
mid- and far-infrared continuum. We compare these relations to similar relations in starburst galaxies in order to 
examine the contribution of AGNs to the bolometric luminosities of their host galaxies. The spectra are available 
CZJ . in electronic form from the authors. 

> 

• 1. Introduction Polycyclic Aromatic Hydrocarbons) as a new tool for test- 

ing unified schemes in Seyfert 1 and Seyfert 2 galaxies, and 
Mid-infrared spectra contain a large number of atomic, ^ ^ &1 (200Qa) haye searched for broad line compo _ 

ionic and molecular lines along with various solid-state nentg q{ hidden hiQ&d ^ (BLRs) Fm & ^ 

and dust features in grains of different sizes. They are , , , i q r< i /onnn\ a 

„ , y complete, recent review see Genzel & Gesarsky (2000). A 

therefore an excellent tool to study the nature of infrared , . e , , n , , . ,. r A ^ TO \rv\j 

J good overview of the global characteristics of MIR AGN 

bright galaxies, such as galaxies with active galactic nu- gpectra (nke shapeS; featureg) relatiye ^ intensities) can 

clei (AGNs). For instance, mid-infrared (MIR) continua fee obtained from the complete 2 5 to 45 ^ spectra of 

can be compared to torus models (e.g. D. Alexander et geveral AQN &nd starburgt prototvpes presented in Sturm 

al. 1999, Siebenmorgen et al. 1997), whereas ice absorp- ^ ^ (200Q); ^ combination with the detailed and com _ 

tion features allow the study of the composition of the hensive line spectra of NGC 1068; NGC 4151 and the 

ISM and the extinction law in AGNs (Chiar et al. 2000, n- ■ , , r , , • T . . , /r, nnfl , \ c . . 

. v , ' Gircmus galaxy published m Lutz et al. (2000b), Sturm et 

Spoon et al. 2002). Clavel et al. (2000) have used aro- ^ (ig9g) and Moorwood et aL (ig96) 
matic infrared bands (AIBs, often identified as 'PAHs' - 

Send offprint requests to: sturm@mpe.mpg.de In this P a P er we concentrate on narrow emission lines. 

* Based on observations with ISO, an ESA project with in- The MIR fine structure emission lines cover a wide range 

struments funded by ESA member states (especially the PI of physical conditions, like excitation (0-300 eV), or criti- 

countrics: France, Germany, the Netherlands and the United cal density (10 2 < n H < 10 7 cm" 3 ). They are significantly 

Kingdom) and with participation of ISAS and NASA. less affected by extinction than UV, optical or near-IR 
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Table 1. The sample 
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Nr. 


Galaxy 


Type 


z 


-- 

°12 


-- 

^25 


°60 


°100 


D ft 


_ 4- + + 


1 


Cen A 


Sy2 + SB 


0.001828 


13.3 


17.3 


162 


314 


4 


0.9 


2 


Circinus 


Sy2 + SB 


0.001498 


18.9 


68.4 


248 


314 


3 


0.7 


3 


M 51 


Sy2 + SB 


0.001544 


11.0 


17.5 


109 


292 


6 


1.3 


4 


M 87 


Sy 


0.004360 


0.2 


0.2 


0.4 


1.0 


17 


0.1 





Mkn 1 (=1NCjU 449 J 


by2 


0.015944 


i n 
1.9 


n n 
0.9 


o r 

2.0 


o n 
2.9 


DO 


9.5 


6 


Mkn 3 


Sy2 


0.013509 


0.7 


2.9 


3.8 


3.4 


54 


6.1 


7 


Mkn 6 


Syl.5 


0.018813 


0.2 


0.7 


1.2 


<1.7 


76 


3.6 


8 


Mkn 335 


Syl.2 


0.025785 


0.4 


0.3 


0.3 


0.5 


104 


4.4 


9 


Mkn 463 


(2x) Sy2 


0.050802 


0.5 


1.6 


2.2 


1.9 


202 


51.0 


1U 


Mkn t>09 


Syl.2 


n no/ion? 
0.034,39/ 


A Q 

0.3 


0. ( 


1.4 


1 r 
1.0 


i on 
139 


i a n 
14.0 


11 


Mkn 573 


Sy2 


0.017259 


0.3 


0.9 


1.2 


1.4 


69 


3.5 


12 


Mkn 938 (=NGC 34) 


Sy2 


0.019784 


0.4 


2.4 


16.6 


17.2 


80 


28.0 


13 


Mkn 1014 


Syl 


0.163099 


0.1 


0.5 


2.2 


2.2 


677 


250.0 


14 


NGC 613 


Sy/Radio + SB 


0.004920 


1.0 


2.5 


22 


50 


20 


3.0 


15 


NGC 1068 


Sy2 


0.003829 


38.0 


86.0 


186 


239 


14 


30.0 


16 


NGC 1275 


Sy2 + SB 


0.017559 


1.1 


3.5 


7.1 


6.9 


71 


16.0 


17 


NGC 1365 


Syl.8 + SB 


0.005457 


3.3 


1.1 


76 


142 


22 


12.0 


18 


NGC 3783 


Syl 


0.009730 


0.8 


2.5 


3.3 


4.9 


34 


2.5 


19 


NGC 4051 


Syl.5 


0.002418 


0.9 


1.6 


7.1 


23.9 


10 


0.4 


20 


IN 4151 


Syl.5 


0.003319 


1.8 


4.7 


6.9 


10.2 


20 


1.7 


21 


NGC 5506 


Syl.9/NLXG 


0.006181 


1.3 


3.6 


8.4 


8.9 


23 


2.0 


22 


NGC 5643 


Sy2 


0.003999 


1.1 


3.7 


19.5 


38.2 


16 


1.8 


23 


NGC 7469 


Syl.2 + SB 


0.016398 


1.3 


5.8 


25.9 


34.9 


65 


29.0 


24 


NOn 7589 


Sv9/NI,XO -1- SR 


005954 


1.4 


6.3 


48 


72.8 


20 


5.0 


25 


PKS 2048-57 


Sy2 


0.011348 


1.1 


3.9 


4.3 


4.2 


45 


6.1 


26 


Tol 0109-383 (=NGC 424) 


Syl.8 


0.011661 


1.1 


1.7 


1.8 


1.8 


47 


3.7 


27 


I Zw 1 


Syl 


0.061086 


0.5 


1.2 


2.2 


2.6 


248 


71.0 


28 


I Zw 92 (=Mkn 477) 


Sy2 + SB 


0.037799 


0.1 


0.5 


1.3 


1.9 


153 


12.0 


29 


3C120 (=11 Zw 14) 


Syl 


0.033010 


0.3 


0.6 


1.3 


2.8 


135 


13.5 



t = IRAS (FSC) fluxes at 12, 25, 60 and 100 fim in Jy 
ft = distance in Mpc, for H =75km/s/Mpc and q =0.5 

t f t = 8-IOOOmhi luminosity in 1O 1O [L ], Lir = 5.6xD^ pc (Sioo+2.58S 6 o+5.16S25 + 13.48Si2), see Sanders &Mirabel (1996) 



lines, and their emissivities depend only weakly on elec- 
tron temperature, in contrast to their optical and UV 
counterparts. Low excitation ionic fine structure lines (ex- 
citation potential < 50 eV) sample HII regions that are 
photoionized by massive stars (Spinoglio & Malkan 1992, 
Voit 1992, Giveon et al. 2002). Some contribution from 
ionizing shocks is possible (e.g. Contini & Viegas 2001). 
Ionic lines from species with excitation potentials up to « 
300 eV sample highly ionized or 'coronal' gas and require 
the very hard radiation fields of AGNs, or fast ionizing 
shocks. Electron densities can be determined from ratios 
like [Ne V] 14/24 fim. Line ratios from different ionization 
stages probe the hardness of the radiation field and the 
ionization parameter (e.g. [Nell], [Nelll], [NeV], [NeVI]), 
but possible contributions from nuclear starburst compo- 
nents have to be taken into account. Finally, hydrogen 
recombination lines and line ratios like [S III]33/im/18/im 
can be used to trace the MIR extinction. 

The Infrared Space Observatory (ISO, Kessler et al. 
1996) has provided access to this wealth of information. In 



the following we present a compilation of ISO-SWS emis- 
sion line spectra and line fluxes of AGNs (Sect. |i|). We 
use this data set to derive electron densities and temper- 
atures i n the narrow line regions (NLRs) of AGNs (Sects. 
|5.l| and 5^), and we compare the line ratios of highly 



ionized lines to the predictions of photoionization mod- 
els to constrain the spectral energy distribution (SED) 
of the central ionizing source and to tes t their consistency 
with AGN unification models (Sect. 5.3). Furthermore, we 
study line profiles and their various degrees of asymmetry 
and differences to the profiles of their optical counterparts 
(Sect. 5^4). In Sect. || we develop a new diagram to iden- 
tify those sources in our sample where star forming regions 
within the SWS apertures contribute to the fluxes of low 
lying fine structure lines ('composite sources'). In order 
to establish tools for the distinction of AGNs and star- 
bursts in distant dusty galaxies we construct mid-IR line 
ratio diagrams which are better suited for dusty objects 
than their classical optical counterparts (e.g. Baldwin et 
al. 1981; Veilleux & Osterbrock 1987) in Sect. |. Finally, 
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Table 2. Observed fine structure line fluxes and upper limits (in 10 -20 W/cm 2 ). 



Galaxy 


[Si IX] 


[MgVIII] 


[Si IX] 


[MglV] 


[Fell] 


[MgVII] 


[MgV] 


[Aril] 


'^rcst,vac(A tm ) 


2.584 


3.028 


3.936 


4.487 


5.340 


5.503 


5.610 


6.985 


E ion (eV) 1 


303.2 


224.9 


303.2 


80.1 


7.9 


186.5 


109.2 


15.8 


Resolution 4 


100 


120 


150 


130 


250 


230 


230 


190 


Aperture 2 


14x20 


14x20 


14x20 


14x20 


14x20 


14x20 


14x20 


14x20 


Cen A 




<0.5 


<0.4 




1.6 


1.5 


0.9 3 


4.9 


Circinus 


1.7 


6.2 


4.9 


1.7 


8.3 


8.0 


6.7 


27.6 


M51 




<0.3 










<1.2 




M87 














<1.0 




Mknl 














<2.0 




Mkn3 


















Mkn6 


<0.6 


<0.5 




<1.5 










Mkn335 


<0.5 


<0.3 


<0.2 






<1.8 






Mkn463 












<0.8 


<0.9 




Mkn509 




0.4 














Mkn573 


<0.8 


0.6 










<1.0 




Mkn938 


<1.4 
















Mknl014 




<1.3 










<1.5 




NGC613 


















NGC1068 


3.0 


11.0 


5.4 


7.4 


5.0 


13.0 


18.0 


13.0 


NGC1275 


<0.7 




<0.2 




<1.0 




<1.3 


0.6 


NGC1365 


<0.5 


<0.3 


<0.3 


<1.5 




<1.8 


<1.5 


14.1 


NGC3783 














<1.0 




NGC4051 


















NGC4151 


0.2 


0.6 


0.4 


0.3 




<1.0 


<1.5 




NGC5506 


<1A 


<0.4 


<0.3 


<1.8 


<2.0 


<2.4 


<1.8 


3.2 


NGC5643 














<2.0 




NGC7469 


<0.8 


<0.2 


<0.4 


<1.0 




<2.3 


<2.5 


6.7 


NGC7582 


<0.5 


<0.3 


<0.3 


<1.2 


<1.5 


<0.9 


<1.5 




PKS2048-57 




0.4 


0.2 






<1.4 


<1.4 




TOL0109-383 




<0.2 


<0.2 






<1.1 






I Zw 1 














<0.6 




I Zw 92 


<0.4 






<1.5 










3C120 


<0.4 


<0.4 




<0.8 











1 = lower ionization potential of the stage leading to the transition 

2 = in arcsec 

3 = uncertain detection, either because of noise (~ 3<r) or offset in wavelength larger than 100 km/s compared to similar lines 
in the same object 

4 = in km/s 



we discuss correlations of mid-IR lines with far-IR fluxes in 
the context of the far-IR energy source and the starburst- 
AGN connection (Sect. ||). The rich spectrum of molecular 
hydrogen (H2) in our data set is discussed in a separate 
paper (Rigopoulou et al. 2002). 

2. Target Selection 

Most of the objects discussed here have been observed as 
part of a guaranteed time ISO project on bright galac- 
tic nuclei (MPEXGAL), expanded by a solicited follow- 
up proposal (ZZCORLIN). They are well studied, nearby 
Seyfert galaxies, selected only on the basis of their bright- 
ness and visibility to ISO. The sample is therefore not 



complete in any sense. From the ISO archive we have 
added 9 targets which have been observed in other pro- 
grammes, but for which fluxes have not been published 
to date. For three of these additional sources relative line 
strengths have been previously published and studied by 
Spinoglio et al. (2000). Table lists the names and clas- 
sifications of all 29 galaxies of our sample. According to 
the classifications given in the NASA/IPAC Extragalactic 
Database (NED) the sample consists of 13 Seyfert Is 
(of various sub-types) and 14 Seyfert 2s. NGC 5506 and 
NGC7582 are classified as NLXGs. No Seyfert type is 
given for M 87 and NGC 613. For 8 additional Seyfert 
galaxies (from the CfA sample) SWS fluxes or flux limits 
of 4 coronal lines can be found in Prieto & Viegas (2000). 
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Table 2. continued 
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v. i ciicL^v y 


[Ne VI] 


[Fe VII] 


[ArTTTl /TMp- VTT1 


[Fe VII] 


[S TV1 


[Ne II] 


[Ne V] 


[Ne III] 


fP TTT1 
L in] 


[Fe TT1 
[re iij 




7.652 


7.815 


8.991 


9.527 


10.551 


12.814 


14.322 


15.555 


17.885 


17.936 


E ion (eV) 1 


126.2 


99.1 


27.6/186.5 


99.1 


34.8 


21.6 


97.1 


41.0 


19.8 


7.9 


TT» 1 j ■ 4 

Resolution 


165 


160 


150 


140 


115 


190 


160 


150 


125 


125 


irpci Li Lll c 


14x20 


14x20 


14x20 


14x20 


14x20 


14x27 


14x27 


14x27 


14x27 


14x27 


Cen A 


2.0 




0.9 




1.4 


22.1 


2.7 


14.1 


0.4 


0.2 


Circinus 


36.3 


1.5 


8.3 


1.1 


12.7 


90.0 


31.7 


33.5 


1.3 3 


<1.0 


M51 


<0.3 




<0.5 




<0.8 


7.0 


<0.2 


2.9 






M87 






<0.4 
















Mknl 






<0.9 
















Mkn3 


4.5 










4.7 


4.6 


12.3 






Mkn6 


<0.5 










0.8 


<0.6 


1.8 3 






Mkn335 


<0.3 










<0.7 


<0.5 


<0.8 






Mkn463 






0.3 






1.3 


1.4 








Mkn509 


0.4 










2.0 










Mkn573 


1.2 










<1.3 


1.8 


2.4 






Mkn938 












5.7 


<0.4 


<1.0 






Mknl014 












0.4 3 


<0.4 








NGC613 


<0.3 










4.0 


<0.8 


<6.8 






IN vJX V,' 4UUO 


110.0 


3.0 


25.0 


4.0 


58.0 


70.0 


97.0 


160.0 




<10.0 


NGC1275 


<0.5 




<0.5 




<0.3 


2.9 


<0.7 


1.0 


0.4 


<0.7 


NGC1365 


2.5 




1.1 




2.6 


40.9 


2.5 


7.7 






A T A 1 /" 1 ,1-7,111 

NGC3783 






0.5 3 








2.1 








NGC4051 












2.5 


<1.0 


3.0 








7.8 




2.2 




11.3 


11.8 


5.5 


20.7 






NGC5506 


2.2 




0.7 3 




5.4 


5.9 


2.6 


5.8 






NGC5643 






0.9 
















NGC7469 


1.1 




0.5 3 




0.9 


22.6 


<1.5 


2.2 






NGC7582 


2.3 




1.0 




1.8 


14.8 


2.2 


6.7 




<0.6 


PKS2048-57 


2.8 




1.0 






2.1 3 


3.6 


5.9 






TOL0109-383 


1.2 










<1.2 


0.8 


<1.6 






I Zw 1 


<0.3 










0.65 


0.27 








I Zw 92 


<0.7 










2.4 


<1.5 


1.6 






3C120 


1.7 










0.9 


1.9 


2.0 







Many active galaxies exhibit intense star formation in 
the (circum-)nuclear region. In these cases some fraction of 
the measured fluxes of low lying fine structure lines (exci- 
tation potential < 50 eV) will be produced by photoioniza- 
tion from stars rather than from the AGN. In Table |l| we 
have indicated (Type='Sy+SB') such composite sources 
where star formation regions within the SWS apertures 
are important. This classification will be discussed in Sect. 

I 

3. Observation and Data Processing 

The data have been obtained with the Short- Wavelength 
Spectrometer (SWS, de Graauw et al. 1996) on board the 
Infrared Space Observatory (ISO), using the SWS02 and 
SWS06 observation modes. These modes provide line spec- 
tra with a resolving power of approximately 1000-2500 
in the range between 2.4 and 45 /im. For Circinus and 
NGC 1068 a few lines have been added from SWS01 mea- 



surements with reduced resolving power. The observations 
were centered on the nuclei. Total observation times per 
source ranged from 10 minutes to 7 hours. Spectral lines 
at different wavelengths are observed with three different 
aperture sizes, varying between 14"x20" and 20" x 33" (see 
Table |). 

We have processed the data using the SWS Interactive 
Analysis (IA) system (Lahuis et al. 1998, Wieprecht et 
al. 1998) and the ISO Spectral Analysis Package ISAP 
(Sturm et al. 1998) together with the associated calibra- 
tion files available at the end of 2000/beginning of 2001. 
Dark current subtraction, scan direction matching, and 
flatfielding have been done interactively, and noisy detec- 
tors have been eliminated. Outliers were clipped, and the 
data of all 12 detectors were averaged for each line, re- 
taining the instrumental resolution. In a few cases some 
wavelength ranges were affected by fringes. In these cases 
we have defringed their averaged spectra using the FFT or 



Table 2. continued 
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Galaxy 


[S III] 


[Fe III] 


[Fel] 


[Ne V] 


[SI] 


[0 IV] 


[Fe II] 


[S III] 


[Sill] 


[Fe II] 


[Nc III] 


Arcst,vac(/^m) 


18.713 


22.925 


24.042 


24.318 


25.249 


25.890 


25.988 


33.418 


34.815 


35.349 


36.014 


E ion (eV) 1 


23.3 


16.2 


0.0 


97.1 


0.0 


54.9 


7.9 


23.3 


8.2 


7.9 


41.0 


Resolution 4 


120 


250 


230 


230 


225 


215 


215 


195 


185 


185 


185 


A 2 

Aperture 


14x27 


14x27 


14x27 


14x27 


14x27 


14x27 


14x27 


20x33 


20x33 


20x33 


20x33 


Cen A 


6.4 


0.3 


<0.2 


2.0 


<0.3 


9.8 


1.2 


22.3 


54.5 




1.4 


Circinus 


35.2 


1.2 


<1.0 


21.8 


<0.5 


67.9 


5.9 


93.2 


151.0 


1.5 


4.0 


M51 


1.0 










1.9 


1.0 


4.6 


9.6 






M87 


<0.6 










<0.2 


0.5 


<1.4 


<1.2 






Mknl 


<0.8 










2.8 


<0.4 


^1 


<"1 3 

\ J- .o 






Mkn3 








3.4 




12.6 




3.0 






<2.5 


Mkn6 








<0.5 




1.6 3 




<1.8 








Mkn335 








<0.3 




1.3 


<0.3 








<0.5 


Mkn463 


<0.8 










4.3 




1.2 


2.0 






Mkn50Q 


<1.5 






<0.6 




1.8 










<2.0 


Mkn573 












7.9 




2.9 








Mkn938 
























Mknl014 


<0.6 










2.4 


<1.0 


<1.6 


<1.3 






AT S~i f~if> 1 O 

NGC613 








<0.7 














<1.0 


noc!1068 

IN vJX V,' 1UUO 


40.0 






70.0 




190.0 


go 3 


55.0 


91.0 


<rl 4 


18.0 


NGC1275 


<1.3 


<0.5 


<0.4 


<0.4 


<0.5 


<0.5 


<0.4 


<1.7 


6.8 




<1.0 


TVT /~i /~1 -I o r> f 

NGC1365 


13.5 


0.8 


<0.7 


3.9 


<0.4 


14.6 


1.8 


36.1 


73.8 




<1.6 


AT /~i no rroo 

NGC3783 


2.0 










3.8 


<0.5 


<2.0 


3.2 






NGC4051 
























NGC4151 


5.4 


<0.2 


<0.4 


5.6 




20.3 


0.4 


8.1 


15.6 




3.5 


NGC5506 


3.0 


<0.4 


<0.5 


2.4 




13.5 


0.6 3 


8.1 


14.2 




1.2 3 


NGC5643 


2.3 










10.3 




5.8 


5.5 






NGC7469 


9.2 


0.9 


<0.4 


0.6 




3.4 


0.6 


10.4 


19.6 




<1.0 


NGC7582 


5.2 


<0.5 


<0.5 


3.1 




11.6 


0.8 


11.3 


21.8 


<1.4 


1.0 3 


PKS2048-57 


2.2 






2.7 




8.5 


0.5 


2.7 


3.5 




<1.2 


TOL0109-383 








0.5 




1.4 


<0.3 








<2.0 


I Zw 1 


<0.5 










<0.6 




<1.0 


<1.8 




<0.4 


I Zw 92 












1.1 




1.5 








3C120 












7.5 


<0.3 


<1.2 









Table 3. Additional fine structure line fluxes and upper limits (in 10 20 W/cm 2 ). 



Galaxy 


[CalV] 


[CaVII] 


[CaV] 


[ArVI] 


[Si VII] 


[Nalll] 


[NaVI] 


[CaV] 


[MgV] 


[Fell] 


[Felll] 


Arcst,vac(/^ni) 

Resolution 4 


3.207 


4.087 


4.158 


4.530 


6.492 


7.318 


8.611 


11.482 


13.521 


24.519 


33.038 


125 


150 


140 


135 


215 


180 


150 


110 


185 


220 


200 


E ion (eV) 1 
Aperture 2 


50.9 


108.8 


67.3 


75.0 


205.1 


47.3 


138.4 


67.3 


109.2 


7.9 


16.2 


14x20 


14x20 


14x20 


14x20 


14x20 


14x20 


14x20 


14x20 


14x27 


14x27 


20x30 


Circinus 


<0.9 


<3.5 


<3.0 


3.1 


<1.0 






<0.8 


<1.5 


<0.5 


<5.0 


NGC1068 


1.3 






15.0 




5.8 


<16.0 











1 = lower ionization potential of the stage leading to the transition 

2 = in arcsec 



iterative sine fitting options of the defringe module within 4. Results 
ISAP. 

The spectra are displayed in Figs. [l7] to |23|. The fine struc- 
ture emission lines are generally resolved, having widths 
(FWHM) in the range between 200 and 600 km/s, i.e. 
line widths typical of narrow line regions of Seyfert galax- 
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Table 4. Observed hydrogen recombination line fluxes 
(in 1(T 20 W/cm 2 ). 



Galaxy 



Br/? Pf 7 



Bra 



Pf (3 



Pi a 



A r < 

Resolution 1 
Aperture 2 
Cen A 
Circinus 
M51 
Mkn335 
Mkn463 
Mknl014 
NGC1068 
NGC1275 
NGC1365 
NGC3783 
NGC4151 
NGC5506 
NGC5643 
NGC7469 
NGC7582 

PKS2048-57 
TOL0109-383 
I Zw 1 



2.625 

145 

14x20 

0.9 

3.2 

<0.8 

<0.3 

<0.4 

<0.3 

4.1 

<0.2 

1.6 

<0.8 

0.5 

0.7 4 

0.4 s 

0.75 6 

0.9 

<0.5 
<0.5 
<0.4 



3.741 

155 

14x20 



<0.4 



4.052 

150 

14x20 

0.8 

5.5 



4.654 

135 

14x20 



6.9 

<2.2 

2.8 

0.9 3 
1.2 



1.6 



<10.0 <3.0 



0.6 



7.460 

180 

14x20 

0.2 

2.1 



1 = in km/s 

2 = in arcseconds 

3 = from Br (3 (case B, see Hummer & Storey 1987) 

4 = from Bra(see note 2) 

5 = from Br 7 (Kowara et al. 1987, see note 2) 

6 = from Br 7 (Goldader et al. 1995, see note 2) 



ies. Measured line fluxes and upper limits are summarized 
in Table § (fine structure lines), Table |] (additional fine 
structure lines for few sources), and Table [| (hydrogen re- 
combination lines). The spectra and fluxes for NGC4151 
and NGC 1068 have already been published in Sturm et 
al. (1999) and Lutz et al. (2000b) respectively. We have 
repeated the fluxes (but not the spectra) here for conve- 
nience. Some of the Circinus data have been published in 
Moorwood et al. (1996), but later on more data have been 
taken. We have re-analyzed the old Circinus data with the 
newer calibration files in order to be consistent with the 
reduction of the additional observations. 



5. Physical conditions in the NLR 

Diagnostics of the physical conditions in the Narrow Line 
Region are of interest given the wide range of possible ex- 
citation mechanisms and interpretations for its emission. 
These include models assuming ionization-bounded pho- 
toionization regions (e.g. T. Alexander et al. 1999, 2000), 
more complex photoionization models invoking a mix- 
ture of ionization-bounded and matter-bounded clouds 
(Binette et al. 1996, 1997) as well as shock excitation, for 
example, through jet-cloud interaction (e.g. Axon et al. 
1998). Traditional optical diagnostics have been applied to 
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Fig. 1. The ratio of [NeV] 14.3/24.3 fxm. Diamonds: 
Seyfert Is. Stars: Seyfert 2s, Plus signs: NLXGs. The 
theoretical low density limit («0.9) and the ratio for 
n e =10 4 cm~ 3 (2.7), for a temperature of 10000K, are 
shown as dashed lines. Error bars for a 20% line flux error 
are indicated. The numbering of galaxies can be found in 
the first column of Table [|. 

large sets of objects but do have shortcomings. The [S II] 
AA6716,6731 density diagnostic applies to low ionization 
species that may not be representative of the higher ex- 
citation regions and may be contaminated by non-AGN 
emission. The [OIII] (A4959 + A5007)/A4363 temperature 
diagnostic includes a weak and blended line that is no- 
toriously difficult to measure. Combining infrared fine- 
structure lines with optical lines provides additional op- 
tions. Given the large observing apertures for mid-IR in- 
struments like ISO-SWS, their use should be restricted 
mainly to high excitation species that are not contami- 
nated by spatially extended starburst emission. 

5.1. Density 

Density-sensitive ratios of bright infrared fine-structure 
lines are observable for a range of ionization potentials. 
Because of possible contamination by starburst emission 
(cf. Lutz et al. 2000b for the longer wavelength lines in 
NGC 1068), we do not discuss [S III] 18.71/33.48/zm and 
[Nelll] 15. 55/36. 01/im. The perhaps most interesting ra- 
tio is [NeV] 14.3/24.3 /*m. These lines are not diluted by 
starburst contributions, since they have a lower ioniza- 
tion potential of 97eV and are undetected in ISO spectra 
of starburst galaxies (Genzel et al. 1998). Their ratio is 
largely insensitive to electron temperature and to extinc- 
tion variations (see T. Alexander et al. 1999, their Fig. 3, 
for the theoretical [Ne V] ratio as function of density) . 

All objects with measurable [NeV] in our sample have 
[NeV] ratios indicating a low density (n e = a few 100 to 



a few 1000 



,-3 



well below the critical densities (« 5 
x 10 4 and 5 x 10 5 cm" 3 ), see Fig. [l]. The average ratio 
is 1.1+0.4 for the (4) Seyfert Is, and 1.3+0.3 for the (6) 
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Fig. 2. Temperature-sensitive ratio of the 3426 A and 
24.3/im lines of [NeV] as a function of electron density 
and temperature. The dashed line indicates the ratio for 
NGC 1068 and the dotted lines its uncertainty range con- 
sidering measurement errors and uncertainty of extinction. 



Seyfert 2s, i.e. NLR densities for the two types agree for 
this fairly small sample, consistent with unification. 

5.2. Temperature 

Temperature diagnostics for photoionized gas invoke a 
combination of transitions that originate in levels sepa- 
rated by at least several 1000K. Combinations of ground 
state infrared fine-structure lines with optical/UV for- 
bidden lines of the same ion provide maximum leverage 
for this task and represent excellent diagnostics provided 
practical requirements can be met: (1) Both IR and op- 
tical/UV lines have to be well observed, with decent S/N 
and line-to-continuum. (2) Aperture corrections have to be 
understood - note the much larger ISO aperture compared 
to commonly used optical apertures. (3) Contamination 
by non-NLR emission has to be low. (4) Extinction cor- 
rections have to be understood. The last requirement is 
not trivial to meet. The extinction towards a NLR will 
not follow a simple foreground screen, as most evidently 
shown by the presence of intra-NLR dust (e.g. Cameron 
et al. 1993), and line profile variations between optical 
and infrared (cf. Lutz et al. 2000b for NGC 1068). Lines 
observed at different wavelengths may therefore not sam- 
ple the same gas volume, even if fluxes are adjusted using 
simple extinction corrections. Hence we ignore combina- 
tions UV/IR (e.g. [OIV], [MgV]) where these problems 
are most severe, and discuss a number of optical/IR di- 
agnostics focussing on the well studied objects NGC 1068 
and Circinus. Of particular interest are highly ionized re- 
gions not well sampled by the optical [O III] or [N II] tem- 
perature diagnostics. 
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Fig. 3. Temperature-sensitive ratio of the 6087A and the 
9.53/im (top) and 7.81/xm (bottom) lines of [FeVII], as a 
function of electron density and temperature. The dashed 
line indicates the ratio for NGC 1068 and the dotted lines 
its uncertainty range considering measurement errors and 
uncertainty of extinction. 



One interesting species is [NeV] (lower ionisation po- 
tential 97eV) which has very bright optical and mid- 
infrared transitions. Figure || shows the ratio predicted 
for the 3426A and 24.3/im lines, using the atomic data of 
Lennon & Burke (1994) and Nussbaumer & Rusca (1979). 
The range of allowed values for NGC 1068 (combining our 
data with optical data of Marconi et al. (1996) and adopt- 
ing Ay=0.71) is also indicated, adding a 20% measure- 
ment error for each line, and a 0.25mag uncertainty in Ay. 
to produce the maximum deviation. These uncertainties 
(particularly the one on extinction) lead to a wide range of 
electron temperatures (~17000-45000K) being consistent 
with the NGC 1068 data for its density of 2000 cm- 1 . 
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Another high excitation species offers the potential of 
reducing the extinction effects. The brightest optical line 
of [FeVll] is found at the longer wavelength of 6087A. 
With a lower ionisation potential of 99eV, [FeVII] sam- 
ples a region similar to [NeV]. The ratio of the optical 
[FeVII] 6087A line to the mid- infrared [FeVII] lines is lit- 
tle sensitive to electron density in the regime determined 
above, and forms a diagnostic of electron temperature in 
the NLR. Figure || shows the ratios from solving the rate 
equation using the atomic data of Berrington et al. (2000, 
and priv. communication]^ . We have observed the mid-IR 
[FeVII] lines in NGC 1068 and Circinus. The extinction 
corrected line ratios for NGC 1068 lead to very different 
electron temperatures: «43000K from 6087A/9.53^m but 
«16000K from 6087A/7.81/mi. This discrepancy is also re- 
flected in the extinction corrected ratio of the 9.53/xm and 
7.81/mi lines which is 1.4, significantly different from the 
~3-3.5 expected for a wide range of conditions from the 
atomic data. For Circinus, the extinction corrected ratio 
of the 9.53/Ltm and 7.81^m lines is similarly low, about 1.1. 
There is no immediate explanation for this inconsistency 
between mid-IR FeVII fluxes that is observed indepen- 
dently in two sources. An explanation by extinction uncer- 
tainties is unlikely, since «10mag of additional visual ex- 
tinction would be needed to selectively weaken the 9.53/^m 
line which is inside the silicate feature. Uncertainties in the 
atomic data might be another possibility but are only par- 
tially supported by observations of the same lines in the 
planetary nebula NGC 7027: Salas et al. (2001) observe a 
ratio of 1.84 for the two mid-IR lines, closer to the value 
of 2.73 expected for the conditions in this nebula. 

We summarize that, while there are indications for 
high (20000-30000K) electron temperatures in the part of 
the narrow line region sampled by Ej on ^100eV species, 
there is considerable uncertainty on extinction corrections 
and atomic data which prevents firm conclusions. 

5.3. Excitation 

Varying relative line strengths of low and high ionizaton 
lines, e.g. of the Neon sequence, seem to indicate a large 
variation in excitation among the galaxies in our sample. 
Line ratios from ions in different stages of excitation along 
with a photoionization model can be used to reconstruct 
the NLR radiation field produced by the central ionizing 
source. Mid-Infrared lines are particularly suited because 
they are little sensitive to extinction and electron tempera- 
ture, and because they span a wide range of ionization po- 
tentials. In detailed studies we have analyzed three nearby 
Seyfert nuclei of our sample with elaborate photoioniza- 
tion models: Circinus (Seyfert 2 plus starburst), NGC4151 
(Seyfert 1.5) and NGC1068 (Seyfert 2), see Moorwood et 
al. (1996), and T. Alexander et al. (1999, 2000). These 
models were able to reconstruct the intrinsic spectral en- 
ergy distribution (SED) of the ionizing source in the ex- 
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Fig. 4. Line ratio diagram of typical NLR tracers. 
Diamonds: Seyfert Is. Stars: Seyfert 2s, Plus signs: 
NLXGs. Numbering of galaxies as in Fig.0. The supernova 
remnant RCW103 (Oliva et al. 1999) is shown as a trian- 
gle. The [Ne VI]/[0 IV] values for three additional sources 
are indicated as horizontal arrows on the right hand y- 
axis (from top to bottom: 3C120, Mkn509, Mkn573). To 
constrain ionization parameter and density we have over- 
plotted a model grid by Spinoglio et al. (2000) for a power 
law ionizing continuum of index a = -1.0: solid lines are 
log U = -2.5, -2, and -1.5, dashed lines are for n = 10 2 , 
10 3 , and 10 4 . 



treme UV, where a 'Big UV Bump' around lOOeV is ex- 
pected as a signature from a hot thin accretion disk. The 
results are consistent with such a 'Big UV Bump', but also 
suggest for some sources of both types (Seyfert 1 and 2) 
the presence of neutral absorbers between the AGN's ex- 
treme ultraviolet emitting source and the NLR. Such (UV) 
absorbers have been suggested independently by studies of 
UV absorption lines (Kriss et al. 1992, Kraemer et al. 1999, 
2001). This detailed photoionization modeling requires a 
large number of mid-IR lines (of good S /N) supplemented 
with UV/optical/near-IR lines. For most of the galaxies 
in our sample the compilation of lines from ISO and the 
literature is not as complete as for the three examples 
given above. Hence, wc do not attempt a similar modeling 
for them. Instead, we try to constrain the NLR excita- 
tion by directly comparing our observations to standard 
photoionization models from the literature. 

The main input parameters for photoionization models 
are the electron density n e and the ionization parameter 
U, i.e. the number of ionizing photons per hydrogen atom 
at the inner face of the ionized cloud. The ratio of lines 
from ions of similar ionization potential but with differ- 
ent critical densities (e.g. the [NeV] 14/24 ratio) are good 
tracers of the electron densities (see Sect. 5.1). Vice versa, 
the ratio of lines with similar critical density but from 
ions of different ionization potential (like [Ne VI] / [O IV] ) 
is sensitive to the ionization parameter. In Fig. ^ we have 
used this to construct a diagram to constrain n e and U. 
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Class I 


Class II 


Class Ilia 


Class Illb 


optical [O III] profile 
MIR [0 IV] profile 
examples 


symmetric 
symmetric, =optical 
3C120 1 , PKS2048 2 


asymmetric 
symmetric 

Mkn509 2 , NGC5643 2 


asymmetric 
asymmetric, =optical 
NGC3783 2 ,NGC4151, 3 
NGC7469 1 , Tol0109 2 


asymmetric 

asym. /non-gauss., ^optical 
NGCK^^knl 1 , Mkn3 3 , 
Mkn463 2 , NGC5506 3 



1 = optical line profile from Vrtilek & Carleton 1985 

2 = optical line profile from Whittle 1985 

3 = optical line profile from Veilleux 1991 



We have chosen the lines of [Ne V] , [Ne VI] and [0 IV] be- 
cause they are not affected by photoionization by stars and 
because they are generally among the brightest high ion- 
ization lines. For comparison we show in Fig. ^|the location 
of the supernova remnant RCW103 (Oliva et al. 1999), as 
an example for a strong shock source, and a photoioniza- 
tion model grid taken from Spinoglio et al. (2000). This 
model grid was computed for a power law with an ionizing 
continuum of index a=-1.0 and for various ionization pa- 
rameters and electron densities. We draw two conclusions 
from this comparison: firstly, all galaxies are consistent 
with these standard photoionization models, with average 
ionization parameters log U between -1.5 and -2 and (as 
seen already in Sect. [O]) average densities between a few 
100 and a few 1000 cm -3 ). Such a simple comparison can 
not, however, distinguish between simple power laws and 
more complex 'Big UV Bump' models. As noted earlier, 
the detailed modeling of single sources required for such 
an analysis is outside the scope of this paper. Secondly, 
within our (small) sample, we do not see significant dif- 
ferences between the AGN sub- types. 

5.4. NLR line profiles 

Structures and velocity fields in the NLR can be stud- 
ied by an analysis of line profiles. Numerous studies of 
the emission line profiles in Seyfert galaxies have shown 
that in many cases the lines exhibit blueward asymme- 
tries and are blueshifted with respect to the galaxies' sys- 
temic velocity. The common interpretation for these pro- 
file asymmetries is that they are caused by differential 
extinction in an outflow or inflow of clouds with a modest 
amount of mixed-in dust. Observations in the infrared can 
obviously test these scenarios, because infrared lines suf- 
fer more than an order of magnitude less extinction than 
in the optical. Hence, they should not or only marginally 
show asymmetries. Sturm et al. (1999) and Lutz et al. 
(2000b) have presented the first such studies of optical-to- 
mid-IR line profile comparisons for NGC 4151 and NGC 
1068. In these two sources the outflow-plus-dust scenario 
seems to be (at least partially) wrong, since the MIR lines 
show asymmetries similar to their optical counterparts. 

Optical [O III] line profile information exists for many 
objects in our data set. These can be used for a com- 




-1000 -500 500 1000 -1000 -500 500 1000 



Fig. 5. Comparisons between optical and MIR line pro- 
files (continuum subtracted and normalized to peak flux 
density), for four different cases as discussed in the text. 
Solid lines show the MIR [O IV] 26/im line, taken from the 
ISO-SWS data set. Dashed lines are optical [OIIIJ5007A 
lines, smoothed to the SWS resolution. References for the 
optical line profiles are given in Table [| The velocity scale 
on the x-axis (in [km/s]) is relative to the systemic veloc- 
ities given in Table [I]. 

parison with the MIR [O IV] line, as described in Sturm 
et al. (1999). We have convolved the optical line profiles 
with the SWS instrumental profile (Gauss profiles with 
FWHM as given in Table ||) in order to smooth these pro- 
files to the resolution of SWS. Four such examples of an 
optical-MIR comparison are shown in Fig. ||. The objects 
in our sample can be grouped roughly into three different 
classes in terms of line profile asymmetries and agreement 
of MIR lines with their optical counterparts, as summa- 
rized in Table ||. In some cases (Class I) the optical line 
profile of [O III] is quite symmetric and of gaussian shape. 
In these cases, not surprisingly, the MIR and optical line 
profiles match each other very well. Another group of ob- 
jects (Class II) shows strong blueward asymmetries in the 
optical, while the MIR lines are rather symmetric. This 
is exactly what scenarios with infall our outflow of NLR 
clouds with mixed-in dust predict. We only found two such 
cases in our sample. The third class of objects has asym- 
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metric optical lines, but MIR lines which are inconsistent 
with these scenarios. This class can be further divided into 
two sub-classes: in Class Ilia the MIR lines are asymmet- 
ric, too, and agree well with the optical lines. This case 
has been studied in more detail in an analysis of NGC 
4151 by Sturm et al. (1999), and can be explained, for 
instance, by a true asymmetry in the distribution of the 
NLR clouds, or, in the case of NGC 4151, by a central, 
optically very thick, but geometrically thin absorber on 
parsec scales. In Class Illb the NIR lines are asymmetric 
(or symmetric but with non-gaussian profiles), but differ- 
ent from their optical counterparts. One member of this 
class is NGC 1068, which has been analyzed by Lutz et 
al. (2000b). This suggests that parts of the NLR are sig- 
nificantly obscured in the optical, but not enough to also 
block the MIR lines. Similar to Class Ilia, the remaining 
MIR profile asymmetries may be either due to an intrinsic 
asymmetry of the NLR, or due to a very high density ob- 
scuring component which is hiding part of the NLR even 
from infrared view. We note that for some of the objects 
with MIR lines of good S/N (Cen A, M 51, MKN 573, 
NGC 1365, and NGC 7582) there is, to our knowledge, no 
(suitable) optical line profile information available in the 
literature. 

Many objects in our sample exhibit differences in the 
profiles of lines with different ionization potential. For in- 
stance, the Circinus galaxy shows symmetric low ioniza- 
tion lines, but asymmetric high ionization lines with the 
typical blueward asymmetries. Vice versa, in NGC 7582 
high ionization lines are very symmetric, while the low 
ionization lines are strongly asymmetric. Correlations, as 
well as anti-correlations, of line asymmetries with criti- 
cal density and/or ionization potential have been claimed 
in many publications. For NGC 7582 a contribution 
to the low ionization lines from a starburst component 
with asymmetric spatial distribution could be an addi- 
tional/alternative solution. For the Circinus galaxy the 
situation is even more complex, since Oliva et al. (1994) 
reported asymmetric low ionization lines. Finally, for some 
objects in our sample all MIR line profiles appear to be 
quite similar. For instance, the profiles of Cen A, M 51, 
MKN 573, NGC 1365 (for which no optical counterparts 
exist) are quite symmetric and of gaussian shape, regard- 
less of ionization potential or critical density. 

It appears that there is no unique answer to what 
causes the line profile asymmetries in Seyfert galaxies. 
Extinction by dust on different spatial scales and with 
varying column densities, sub-structure and true asymme- 
tries in the spatial distribution of NLR clouds may all play 
a role with varying degrees of importance. For some ob- 
jects the assignment to a certain group is not unique and 
depends also on S/N. We use this classification purely to 
obtain an overview of this complex issue. We refrain from 
a more detailed analysis of line profiles in this paper be- 
cause it requires a careful study of the effects of aperture 
differences (spatial resolutions) between the optical and 
MIR observations, and, in many cases, a larger data set 
of high ionization lines with good S/N. 
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Fig. 6. The ratio [NeVI]/[OIV] vs. jNe VI] /[Nell]. Same 
symbols and numbering as in Fig. |lL A typical error bar 
(according to a 20% individual line flux error) is shown in 
the upper left corner. An empirical dividing line is drawn 
as dashed line. Composite sources according to their PAH 
spectra are encircled. For galaxies #26, 28 and 29 no ISO 
PAH spectra exist. 

6. Composite sources 

As discussed in Sects. [I] and || for many active galax- 
ies an additional contribution from star forming regions 
within the SWS apertures is to be expected. In these 
cases ('composite sources') some fraction of the measured 
fluxes of low lying fine structure lines (excitation poten- 
tial < 50 eV) will be produced by photoionization from 
stars rather than from the AGN. In the previous sections 
we have therefore concentrated on high excitation lines 
([OIV], [NeV], [NeVI], [FeVII]) with little or no contam- 
ination from possible starburst components. However, for 
our considerations in the following sections low ionization 
lines like [Nell] will be important. Hence, an independent 
identification of composite sources is desirable. 

One such independent indicator is the strength of PAH 
features in the 5 to 12 /jm range, which are known to 
be strong in star forming regions, but more or less ab- 
sent in the spectra of pure AGNs. Genzel et al. (1998), 
Rigopoulou et al. (1999), and Tran et al. (2001) have 
analysed a large number of low resolution ISOPHOT-S, 
ISOCAM-CVF and SWS spectra and used the line-to- 
continuum (L/C) ratio of the 7.7/im PAH feature to dis- 
tinguish between starburst and AGN dominated sources. 
Following their method (and using the ISO archive as 
well as the spectra published in Rigopoulou et al. 1999 
and Spoon et al. 2002) we classify Cen A, Circinus, M 
51, NGC 613, NGC 1365, NGC 7469, and NGC 7582, 
i.e. those sources with a 7.7 /im L/C ratio greater than 
1, as composite sources. As a consistency check, and for 
the cases where no PAH spectrum is available (Tol 0109, 
I Zw 92, and 3C 120), we have plotted in Fig. | the po- 
sition of the targets in a diagram of [NeVI] /[OIV] ver- 
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Fig. 7. The ratio [OIV]25.9Aim / [NeII]12.8/xm. Same 
symbols and numbering as in Fig. [I]. Composite objects 
are encircled. The starburst prototype M 82 is shown as a 
triangle. Denoted on the right hand y-axis is a simple lin- 
ear 'mixing' model of AGN contribution to the bolometric 
luminosity, in percent (see text). The [O IV] /[Nell] ratios 
for 100% AGN and 0% AGN (=100% starburst) are as- 
sumed to be 2.7 and 0.012 respectively. 
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Fig. 8. The ratio [NeV]14.3^m / [NeII]12.8/xm. Same 
symbols and numbering as in Fig. EL Composite objects 
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ple linear 'mixing' model of AGN contribution to the bolo- 
metric luminosity, in percent (see text). The [Ne V]/[NeII] 
ratios for 100% AGN and 0% AGN (=100% starburst) are 
assumed to be =1.1 and respectively. 



sus [Ne VI] /[Nell]. As discussed in Sect. 5.3, the ratio 
[NeVI]/[OIV] traces the AGN excitation without con- 
tamination from a possible star formation component. 
[Ne VI] / [Ne II] on the other hand, which is also an exci- 
tation indicator, can be influenced by contributions from 
star forming regions to the [Nell] line. In such a diagram 
one expects the pure AGNs to lie roughly along a diagonal 
(both ratios increase with increasing hardness of the ra- 
diation). Composite sources, however, will have relatively 
stronger [Ne II] lines and should lie in a region further to 
the left than the pure AGN sources. In fact, all sources 
with known substantial starburst contributions, i.e. with 
strong PAH features, fall onto the left part of the dia- 
gram (Cen A, Circinus, M51, NGC1365, NGC 7469, and 
NGC7582, marked with an enclosing circle in Fig. ||). Pure 
AGNs (like NGC4151 and NGC1068) are located in the 
right hand part. In Fig. [}]we have indicated an empirical 
dividing line (dashed). For simplicity we have chosen the 
line where [O IV] equals [Ne II] . According to this diagram 
one additional source without (good S /N) PAH data, I Zw 
92 (#28), has to be classified as composite source. This 
classification is confirmed by independent UV to near-IR 
observations (Heckman et al. 1997, Gonzalez Delgardo et 
al. 2001). We consider Mkn 509 (#10) a borderline case. 
Another source, NGC1275, can not be shown at its exact 
position in the diagram, but according to the upper limits 



for [O IV] and [Ne VI] (and the measured value for [Ne II] ) 
the position must be in the composite source region. The 
galaxies that we classify as composite sources according 
to this scheme are indicated in Table | as 'Sy + SB'. The 
location of composite sources (as indicated by strong PAH 
features) in the left part of Fig. || suggests that a small 
ratio of [NeVI] and [Nell] (< 0.1) is an independent em- 
pirical indicator for a composite object. 

7. New mid-1 R diagnostic diagrams 

Deep surveys carried out by future infrared missions (such 
as SIRTF, SOFIA, ASTRO-F, or Herschel) will sample in- 
frared bright galaxies over a wide range of redshifts and 
luminosities. Quantitative spectroscopy of mid- infrared 
emission lines will be an important diagnostic tool for de- 
termining the detailed properties of distant, dusty galax- 
ies, the source of the extragalactic background, and the 
origin of nuclear activity in galaxies. Many of the funda- 
mental questions of galaxy formation and evolution de- 
pend substantially on the fraction of the total energy out- 
put of distant sources that is produced by star forma- 
tion rather than AGN activity. A large energy contribu- 
tion from hidden AGNs would complicate the deduction 
of the star formation history of the Universe from galaxy 
luminosity functions. In the past, optical line ratios have 
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been used as a tool for distinguishing between the differ- 
ent possible excitation mechanisms and energy sources in 
galaxies (e.g Baldwin et al. 1981; Veilleux & Osterbrock 
1987). This works very well for objects with dust extinc- 
tion less than Ay^ 5, but becomes unreliable for heavily 
obscured objects (see, e.g., Veilleux et al. 1995; Veilleux 
et al. 1999). Hence, mid-infrared analogues of the classical 
optical diagrams are highly desirable. Genzel et al. (1998) 
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Fig. 11. Another example of MIR diagnostic diagrams - 
here: [OIV] 26^m/Br/3 versus [Fell] 26^m/Br/3. Symbols 
and numbering as in Fig. g. 

presented a first empirical version of such a mid-IR tool, 
a diagram with the flux ratio of a high to a low excitation 
line ([OIV] /[Nell], or [Ne V] / [Ne II] ) on one axis together 
with the AIB (PAH) strength on the other axis. Their find- 
ing that the [OIV] /[Nell] ratio is much higher in AGNs 
than in starbursts can now be confirmed on a broader sta- 
tistical basis. In starburst galaxies this ratio reaches values 
of a few times 0.01 at most (Lutz et al. 1998), whereas in 
the AGNs of our sample it does not drop below 0.1 (1.0 
for pure AGNs), as shown in Fig. |^. For comparison we 
show in Fig. || the ratio of [Ne V]14/im/[NeII]. Those ac- 
tive galaxies in our sample with the lowest ratios in Figs. [?] 
and H are the ones we have identified earlier as composite 
sources (encircled in the figures), i.e. AGNs with a signifi- 
cant contribution from star formation to their spectra. On 
the right hand y-axes of Figs. [| and || we have denoted a 
simple linear 'mixing' model to read off the AGN contri- 
bution to the bolometric luminosity. This model will be 
further discussed in Sect. |[ 

Perhaps the best optical discriminator between pho- 
toionization by power-law spectra and by OB stars is the 
diagram in Fig. 3 of Veilleux & Osterbrock (1987). It em- 
ploys a moderately ionized line ([O III]) on one axis and a 
line of [OI], normalized to hydrogen recombination lines. 
Our data set allows the construction of MIR analogues of 
this 'VO diagram' which are much less prone to extinction. 
We have selected [SIV]10.5/im instead of [OIII], which 
has an almost identical ionization potential, [Si II]34/im in- 
stead of [O I], i.e. a line with an ionization potential below 
13.6 eV, and Br/? for normalization. For comparison we 
have used the ISO-SWS starburst data set of Verma et al. 
(in prep.). Starbursts and AGNs are clearly separated in 
this diagram (Fig. ||), not only because [S IV] is stronger in 
AGNs than in normal starbursts, but also because [Sill] is 
(on average) stronger, due to the higher importance of par- 
tially ionized zones in AGNs. Low metallicity starbursts, 
like IIZw40 and NGC5253, have much harder radiation 
fields than normal starburst galaxies. In these cases [S IV] 
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can be as strong as in AGNs, but [Sill] is much weaker. 
The overall appearance of Fig. ^ is thus very similar to the 
optical VO diagram. In contrast to its optical counterpart 
our mid-IR version also includes Seyfert 1 galaxies, since 
we did not detect any broad line components in the Br /3 
lines. For a complete comparison of the optical and mid- 
IR diagrams our data set is missing galaxies of the LINER 
type. Their position in the mid-IR diagram is expected in 
the lower right corner, but remains to be tested in future 
infrared missions. 

Several other versions of mid-IR diagnostic diagrams 
(with lines at different wavelengths) are similarly well 
suited for a distinction of excitation mechanisms. This al- 
lows to cover different redshift ranges or to adjust the 
method to the wavelength coverage of different detectors. 
For instance, it is possible to replace [SIV]10.5/im with 
[OIV]26^m, and [SiII]34/xm with [FeII]26/im (Figs. |o| 
and y_lj). Similar diagnostic diagrams, based on theoret- 
ical modelling, with different sets of (mostly weaker) lines 
have been proposed by Spinoglio & Malkan (1992) and 
Voit (1992). In contrast to the optical VO diagrams, the 
mid-IR versions can be applied to dusty systems with 
much higher extinctions, such as Ultraluminous Infrared 
Galaxies (ULIRGs). ULIRGs are believed to be local ana- 
logues of those distant dusty galaxies, which are (and will 
be) found in great numbers in deep infrared galaxy sur- 
veys. As discussed above, the identification of their energy 
source will be a major task for future infrared missions. 
In Fig. [n] we show the position of the well known ULIRG 
Arp 220 (Sturm et al. 1996). Its position in the diagram 
is consistent with predominant powering by intense star 
formation, which is in accordance with earlier studies (e.g. 
Sturm et al. 1996, Genzel et al. 1998). 

8. The source of the far-infrared luminosity 

In the previous section we have presented tools to detect 
AGNs in dusty galaxies. The presence of an AGN, how- 
ever, does not automatically imply that the AGN activity 
is the dominant energy source in these galaxies. A sig- 
nificant fraction of the bolometric luminosity of Seyfert 
galaxies is emitted in the mid and far infrared by ther- 
mal (dust) emission. For a quantitative understanding of 
the contribution from AGNs to the total energy output of 
distant galaxies it is therefore essential to know whether 
the source powering the IR emission is (re-)radiated light 
from the AGN or from distributed star formation in the 
host galaxy. 

In the following we will compare luminosities in dif- 
ferent tracers for Seyferts and Starbursts in an attempt 
to shed light on these issues. We have selected [O IV] (as 
a high excitation line), [Nell] (as a low excitation line), 
and the mid- and far-infrared continua. In all cases, we 
will keep track of the composite AGNs already identified 
in Sect. ^. 

The mid-IR (A < 40/im) emission of AGNs is domi- 
nated by dust continuum re-radiating energy from black 
hole accretion. Orientation effects in a unified scheme 



cause variations in its level between Seyfert types 1 and 
2, and extended star formation can contribute, more visi- 
bly in Sy2s with their lower AGN continuum (Clavel et 
al. 2000). In contrast, the origin of the 40-120//m far- 
infrared emission of AGNs is considerably uncertain. Both 
the AGN proper and star formation have been advocated 
as energy source in the literature (cf. the review by Genzel 
& Cesarsky 2000 and references therein). For starbursts, 
the mid-IR emission will be dominated by the aromatic 
emission features and transiently heated small grains. The 
far-IR emission originates in large dust grains, and all 
are heated by newly formed or older stars (e.g. Desert, 
Boulanger & Puget. 1990). 

[OIV] originates purely from the NLR in AGNs. In a 
unified scheme, the NLR line luminosity should be orien- 
tation independent and a good tracer of AGN power, in 
particular when using an extinction insensitive and mod- 
est excitation line like [O IV] . The much fainter starburst 
[O IV] emission relates to shocks or some very hot stars 
(Lutz et al. 1998). [Nell] is a pure and fairly good tracer 
of the hot star emission in starbursts. In AGNs the [Nell] 
from the NLR is more easily contaminated by starburst 
emission than the higher excitation [O IV] line. 

We have computed mid-IR (12+25^im) and far- 
IR (60+lOO^m) band luminosities from the respec- 
tive IRAS fluxes following Sanders & Mirabel (1996), 
using Fmir = 1.26x(5.16S 25 +13.48Si 2 )xl0- 18 W/cm 2 , 
and F Fm = 1.26x(Sioo+2.58S 60 )xlO- 18 W/cm 2 , where 
S is the flux density in Jy (from IRAS FSC), and 
L=4ttD 2 F = 3.132xl0 4 xD 2 4pc F L Q . Our results are sum- 
marized in Figs. ^2[]l5| , showing correlations of [O IV] and 
[Nell] line luminosities with MIR and FIR continuum 
luminosities, for pure Seyferts, composite Seyferts (both 
from this paper) and starbursts (from Verma et al. in 
prep.). The object classes are discriminated by plot sym- 
bols. For each object class separately, good correlations 
are observed (see also Prieto et al. 2002 showing a simi- 
lar analysis for a much smaller Seyfert sample and using 
ISOPHOT fluxes). The scatter is partly intrinsic but will 
also include a component due to aperture differences be- 
tween SWS lines and IRAS continua. This is irrelevant for 
the AGN NLRs which are smaller than the SWS aperture 
but may affect some of the more extended starburst emis- 
sion. The slopes are near 1 suggesting that the spectral 
properties of each object class do not change much over the 
overall luminosity range probed by our sample. Comparing 
different object classes, both Fig. [l2] and Fig. |l4| clearly 
reflect the much stronger [O IV] emission of Seyferts com- 
pared to starbursts (Genzel et al. 1998, Lutz et al. 1998). 
The composite Seyferts tend to fall between pure Seyferts 
and starbursts, again confirming our classification based 
on PAH emission, and consistent with the fine structure 
line ratios (Sects. As seen in Fig. [l^ and already 

noted by Genzel et al. (1998), the ratio of [Nell] and FIR 
luminosity is very similar in Starbursts and AGNs. The 
difference in [Nell] /MIR between Seyferts and starbursts 
is somewhat larger (Fig. |l3|), due to the moderately MIR- 
brighter AGN SEDs (e.g., de Grijp et al. 1985, 1987). The 
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Fig. 12. [OIV] luminosity vs. MIR luminosity for AGNs Fig. 13. [Nell] luminosity vs. MIR luminosity for AGNs 
(diamonds) and starburst galaxies (stars). Composite (diamonds) and starburst galaxies (stars). Composite 



sources (see text) are indicated as triangles. 



sources (see text) are indicated as triangles. 



small variations in the ratio of [Nell] and infrared lumi- 
nosities make this line a useful reference line for compar- 
ison of line ratios between different types of galaxies (see 
below). 

A possible test for the origin of the AGN far-infrared 
emission is to compare the quality of the correlations 
between [OIV] (assumed to trace AGN luminosity) and 
MIR/FIR. If the AGN MIR emission were AGN heated 
but not the FIR emission, one would expect less scatter 
in the ratio [OIV] /MIR than in [OIV] /FIR. Figure [jj] 
indeed suggests the correlation with mid-infrared emis- 
sion to be better, but the effect is not strong (see also 
the similar conclusion of Prieto et al. 2002). Partly, the 
small magnitude of the effect may be simply due to the 
scatter in the 'reference' ratio [OIV] /MIR. Even if [OIV] 
were a perfect measure of AGN power (which it clearly is 
not), the orientation effects and larger scale contributions 
to the mid-IR AGN emission (Clavel et al. 2000) would 
cause some scatter in the left panel of Fig. 16. These re- 



sults are consistent with a scenario where the mid-infrared 
luminosity in AGNs is dominated by the active nucleus, 
whereas the far-infrared luminosity is mainly produced 
by star formation. The considerable correlation between 
[OIV] and far-infrared luminosities in Seyferts (Fig. |l4|) 
would then simply reflect a starburst- AGN connection in 
which starburst power and AGN power are linked. There 



is plenty of evidence reported in the literature for such a 
connection. For instance there could be a feedback mech- 
anism that controls (circum-) nuclear star formation and 
Black Hole accretion rate. A more trivial link could sim- 
ply be the common gas consumption habits of starbursts 
and AGNs (see e.g. Cid Fernandes et al., 2001). On the 
other hand, we can not exclude other scenarios in which, 
for example, the FIR luminosity of AGNs is powered by 
the active nuclei. The relation of the coronal line fluxes 
to the mid/far infrared continuum fluxes does not allow 
us directly to quantify the AGN contribution to the total 
luminosity of their host galaxies with sufficient accuracy. 
However, the constancy of the [Ne II]-to-FIR ratio can be 
used to construct models for the AGN/starburst contribu- 
tion to L Boi from the ratios of MIR lines to [Ne II] . In Figs. 
[?] and U (right hand y-axes) we show a simple linear mixing 
model, using the respective median ratios of [O IV] /[Nell] 
and [NeV]/[NeII] for the pure AGNs and the pure star- 
bursts in our sample (see figure caption). Those sources 
that we have identified as composite sources by means of 
their PAH spectrum (or their position in Fig. ^) lie in the 
10 to 50 percent domain of these mixing curves. When ap- 
plying such simple models to individual sources one should 
keep in mind that there is some scatter in the underlying 
relations, and the derived contributions in percent should 
not be taken too literally. For statistical studies, however, 
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these models provide a valuable means of quantifying the 
AGN contribution to the luminosity of dusty, composite 
objects. 

9. Summary 

We have presented an inventory of mid-infrared (MIR) 
lines detected in medium resolution (R~1500) ISO-SWS 
2.4-45/xm spectra of a sample of 29 galaxies with active 
nuclei. This data set is rich in fine structure emission 
lines tracing the narrow line regions and (circum-) nuclear 
star formation regions, and it provides a coherent spec- 
troscopic reference for future extragalactic studies in the 
mid-infrared. From the [NeV] 14/24 /im line ratios we 
conclude that all targets in our sample with measurable 
[Ne V] emission have relatively low NLR electron densities 
(a few 100 to a few 1000 cm -3 ). Electron temperatures are 
difficult to derive. There are indications for high (20000- 
30000K) electron temperatures in the part of the narrow 
line region sampled by Ei on ^lOOeV species, but there 
is considerable uncertainty on extinction corrections and 
atomic data which prevents firm conclusions. 

A comparison of the observations to photoionization 
models showed that all objects are consistent with stan- 
dard photoionization models, i.e. the spectral energy dis- 
tribution of their central ionizing source can be explained 
by a normal power law. These simple comparisons can not, 



however, distinguish between simple power laws and more 
complex models of a 'Big UV Bump' and possibly absorp- 
tion. The detailed modeling required for such an analysis 
is outside the scope of this paper. Within the statistical 
limitations of our sample we do not see any significant 
differences of NLR densities or ionization between the dif- 
ferent AGN types. 

The MIR lines in our data set show various degrees 
of profile asymmetries and different levels of similarities 
to their optical counterparts. It appears that the cause 
for the line profile asymmetries in Seyfert galaxies is a 
complex mixture of extinction by dust on different spatial 
scales, sub-structure and true asymmetries in the spatial 
distribution of NLR clouds. 

We have used the line ratios of [NeVI]/[OIV] and 
[Nc VI] /[Nell] as well as the strength of the PAH fea- 
tures as independent indicators for the identification of 
composite sources, i.e. sources in which star-forming re- 
gions within the aperture contribute to the fluxes of low 
lying fine structure lines. As a new tool for determining 
the properties of dusty galaxies we provide mid-IR line ra- 
tio diagrams, like [OIV]25.9/mi/[NeII]12.8/im, which can 
be used to distinguish between emission excited by active 
nuclei and emission from (circum-nuclear) star forming 
regions. In addition, we present two-dimensional diagnos- 
tic diagrams that arc fully analogous to classical optical 
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Fig. 16. [OIV] luminosity vs. MIR luminosity (left) and 
vs. FIR luminosity (right) for pure AGNs (diamonds) and 
composite sources (triangles). The ±1ct range around the 
average luminosity ratios (computed without composite 
sources) is marked with dashed lines. 



diagnostic diagrams, but better suited for objects with 
high extinction. Several combinations of lines like [O IV] , 
[SIV], [Sill], or [Fell], normalized to hydrogen recombi- 
nation lines, can be used, allowing to adapt to redshift 
range or wavelength coverage of different detectors. 

The high ionization lines like [OIV], [NeV], or [NeVI] 
in Seyfert galaxies are well correlated with both the mid- 
IR and far-IR continuum radiation, the latter correlation 
showing a somewhat larger scatter than the former. Low 
ionization lines, which might be significantly affected by 
star formation, show similar correlations. In particular the 
[NeII]-to-FIR correlation is very similar in both Seyfert 
and starburst galaxies. On the one hand, this makes it 
difficult to quantify the dominant contributor to the bolo- 
metric luminosities of Seyfert host galaxies directly from 
such ratios of MIR lines to FIR continuum. On the other 
hand, it makes [Nell] a valuable reference line for com- 
parisons of line ratios between different sources. We have 
used this to construct simple linear mixing models of AGN 
contribution to the bolomctric luminosities of composite 
sources for the ratios [OIV] /[Nell] and [Ne V] / [Ne II] . 
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Fig. 17. ISO-SWS spectra of Centaurus A and Circinus. 
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Fig. 18. Further SWS spectra. Same units as in Fig. [l7| The range for the combined [O IV] /[Fell] spectra is 
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Fig. 20. Further SWS spectra. Same units as in Fig. [It]. 
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Fig. 21. Further SWS spectra. Same units as in Fig. [I?]. The range for the combined [O IV] /[Fell] spectrum is 
700,+1700]km/s around [O IV]. 
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Fig. 22. Further SWS spectra. Same units as in Fig. [It]. 
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Fig. 23. Further SWS spectra. Same units as in Fig. [Tt]. 
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